Rapid structural change in low-lying states of neutron-rich Sr and Zr isotopes 
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The rapid structural change in low-lying collective excitation states of neutron-rich Sr and Zr iso- 
topes is studied by solving a five-dimensional collective Hamiltonian with parameters determined by 
both relativistic mean-field and non-relativistic Skyrme-Hartree-Fock calculations using the PC-PK1 
and SLy4 forces respectively. Pair correlations are treated in BCS method with either a separable 
pairing force or a density-dependent zero-range force. The isotope shifts, excitation energies, elec- 
tric monopole and quadrupole transition strengths are calculated and compared with corresponding 
experimental data. The calculated results with both the PC-PK1 and SLy4 forces exhibit a picture 
of spherical-oblate-prolate shape transition in neutron-rich Sr and Zr isotopes. Compared with the 
experimental data, the PC-PK1 (or SLy4) force predicts a more moderate (or dramatic) change in 
most of the collective properties around N = 60. The underlying microscopic mechanism responsible 
for the rapid transition is discussed. 

PACS numbers: 21.10.-k, 21.10.Re, 21.60.Jz, 27.60.+J 



I. INTRODUCTION 

The spectroscopy of nuclear low-lying states provides 
rich information about the interplay of nuclear collec- 
tivity and single-particle structure. For neutron-rich Sr 
and Zr isotopes, lots of spectroscopic data have been 
measured [1-11]. The abrupt changing in the lifetimes 
and excitation energies of 2+ states, two-neutron sepa- 
ration energies and rms charge radii indicates a sudden 
onset of large quadrupole deformation at neutron num- 
ber TV = 60. In the past decades, various models have 
been employed to study such a dramatic transition in the 
low-energy structure of these nuclei [12-25], the descrip- 
tion of which has been a challenge in theoretical nuclear 
physics. 

Two main mechanisms based on shell models and 
moan-field approaches were proposed to explain the sud- 
den onset of large nuclear collectivity at TV = 60, namely, 
a strong isoscalar proton-neutron interaction between 
particles occupying the 59/2-57/2 spin-orbit partners [12], 
and the occupation of \ow-K components of the hxi/2 
neutron orbit [26, 27] respectively. Besides these two 
factors, the simultaneous polarization of 2pi/%, 2p 3 / 2 , and 
I/5/2 proton orbits as one goes from 96 Sr to 98 Sr and from 
98 Zr to 100 Zr was pointed out to be the the major factor 
in the more recent projected shell model study [28]. How- 
ever, in most of the previous studies for nuclear low-lying 
states with, for instance the (projected) shell model, ei- 
ther a phcnomcnological effective interaction or effective 
charges for neutrons and protons defined within a specific 
valence space was introduced. In particular, it has been 
shown recently that the large-scale shell model calcula- 
tions with a more extended model space and carefully 
adjusted effective interaction are able to reproduce the 
low-spin spectroscopic data of 90_98 Zr, but fail to give 



large deformation properties for 100 Zr [29]. The mecha- 
nism responsible for the rapid transition in the structure 
of low-lying states in Sr and Zr isotopes around N = 60 
requires further investigation. 

Nuclear energy density functional theory (DFT) is 
nowadays one of the most important microscopic ap- 
proaches for large-scale nuclear structure calculations in 
medium and heavy nuclei. The main ingredient of DFT 
is the energy density functional (EDF) that depends on 
densities and currents representing distributions of nucle- 
onic matter, spins, momentum, and kinetic energy and 
their derivatives. In the past decades, a lot of efforts 
have been devoted into finding out an EDF with reli- 
able and good predictions by optimizing about 10 uni- 
versal parameters to basic properties of nuclear matter 
and some selected nuclei. At present, there are three 
types of most successful EDFs, i.e., the non-relativistic 
Skyrme and Gogny forces and the effective relativistic 
Lagrangian, being employed extensively in the descrip- 
tion of nuclear structure properties [30]. 

The nuclear DFT of single-reference state (SR-DFT) 
with constraint on quadrupole moments, in the context 
of Hartrcc-Fock (HF) or Hartrcc-Fock-Bogoliubov (HFB) 
methods, has already been adopted to study the evo- 
lution of deformation energy surfaces in /3-7 plane (see 
the recent work [31]) for the neutron- rich Sr and Zr iso- 
topes [27, 32, 33]. Even though, all of these studies have 
indeed shown the increasing of deformations based on the 
shift of global minimum. However, most of these stud- 
ies are focused on the properties of nuclear mean-field 
ground states. In particular, a competing weakly oblate 
deformed minima coexisting with a large prolate one has 
been found in nuclei around N = 60, in which case the 
dynamic correlation effects from quadrupole fluctuation 
are expected to be significant. Moreover, the evolution 
behavior of low- lying excited states is sensitive to the bal- 
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ance between the oblate and prolate minima. Therefore, 
it is necessary to extend these studies for the low-lying 
excited states in Sr and Zr isotopes. 

In the context of generator-coordinate method (GCM), 
the framework of SR-DFT has been extended for study- 
ing the nuclear low-lying states. In Ref. [34], the collec- 
tive quadrupole and octupole excitations in Zr isotopes 
were studied using a basis generated by the HF+BCS 
calculations with the SkM* effective interaction. How- 
ever, the low- lying states have not been described quite 
well, probably due to the missing of effects from tri- 
axiality and angular momentum projections. In recent 
years, the GCM has been extended much further by im- 
plementing the exact one-dimensional [35-37] or three- 
dimensional [38-41] angular momentum projection or 
together with particle number projection before or af- 
ter variation for the mean-field states in the modern 
EDF calculations. The dynamic correlation effects re- 
lated to the symmetry restoration and quadrupole fluc- 
tuation (along both /3 and 7 directions) around the mean- 
field minimum are included naturally without introduc- 
ing any additional parameters. However, the application 
of these methods with triaxiality for systematic study 
is still much time-consuming. Up to now, such kind of 
study is mostly restricted to light nuclei [42, 43] and some 
specific medium heavy nuclei [44, 45]. 

As a Gaussian overlap approximation of GCM, the col- 
lective Hamiltonian with parameters determined by self- 
consistent mean-field calculations is much simple in nu- 
merical calculations, and has achieved great success in de- 
scription of nuclear low-lying states [46-50] and impurity 
effect of A hyperon in nuclear collective excitation [51]. 
In particular, a systematic study of low-lying states for 
a large set of even-even nuclei has been carried out with 
the Gogny D1S force mapped collective Hamiltonian and 
good overall agreement with the low-lying spectroscopic 
data has been achieved [52]. However, some fine struc- 
tures along the isotonic or isotopic chains are still not 
reproduced satisfactorily. For Sr and Zr isotopes, the 
evolution of isotope shifts and the change in the proper- 
ties of low-lying states are found to be more moderate in 
comparison with the data. Therefore, in this work, we 
will examine the evolution of low-lying states obtained 
from the calculations with two other popular EDFs, i.e., 
the non-relativistic Skyrme force and the effective rela- 
tivistic Lagrangian. The difference in the results of these 
two calculations will be emphasized. 

The paper is organized as follows. In Sec. II we will 
introduce briefly our method used to study the low-lying 
states in neutron-rich Sr and Zr isotopes. The results 
and discussions will be presented in Sec. III. A summary 
is made in Sec. IV. 



II. THE METHOD 

The quantized five-dimensional collective Hamilto- 
nian (5DCH) that describes the nuclear excitations of 



quadrupole vibrations, rotations, and their couplings can 
be written in the form [46, 47, 49] 



H — T vi b + T n 



Vcoll 



(1) 



where V^ ll is the collective potential. The vibrational 
kinetic energy reads, 
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and rotational kinetic energy, 
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with Jk denoting the components of the angular momen- 
tum in the body-fixed frame of a nucleus. It is noted 
that the mass parameters Bpp, Bp 1 ^ -B 77 , as well as the 
moments of inertia Ik, depend on the quadrupole defor- 
mation variables P and 7, 

I k = 4S fc /3 2 sin 2 (7- 2fc7r/3), fc = l,2,3. (4) 

Two additional quantities that appear in the expres- 
sion for the vibrational energy: r = B1B2B3, and 
w = BppB^j — Bp^, determine the volume element in the 
collective space. The corresponding eigenvalue problem 
is solved using an expansion of cigenfunctions in terms 
of a complete set of basis functions that depend on the 
deformation variables P and 7, and the Euler angles (j), 6 
and ip [53]. 

The dynamics of the 5DCH is governed by the seven 
functions of the intrinsic deformations P and 7: the col- 
lective potential V co \\, the three mass parameters: Bpp, 
Bpy, -B-y-,, and the three moments of inertia Ik- These 
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functions are determined by the relativistic mean-field 
(RMF)+BCS calculations using the PC-PK1 force [54] 
for the particle-hole (ph) and the separable pairing force 
(adjusted to reproduce the pairing properties of the 
Gogny force D1S in nuclear matter) [55, 56] for the 
particle-particle (pp) channels or by the Skyrme-Hartree- 
Fock (SHF)+BCS calculations using the SLy4 force [57] 
for the ph channel and a density-dependent (5-force in the 
pp channel, 



V(v 1 ,r 2 ) = V pp 
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with a strength of Vg p = -1000 MeV fm 3 and p Q = 0.16 
fm~ 3 for both neutrons and protons and with a soft cutoff 
at 5 MeV above and below the Fermi energy as defined 
in Ref. [58] . A constraint on the deformation parameters 
of both P, ranging from to 0.8 (A/3 = 0.05) and 7, 
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FIG. 1: (Color online) Moment of inertia along i-direction 1%, mass parameters Bpp, B 77 and squared proton radii in 5DCH 
for (left panel) 96 > 98 > 100 g r anc j (right panel) 98 > 10 °. 102 Zr as functions of axial deformation parameter /3 determined from the 
mean-field calculations with both the SLy4 and PC-PK1 forces. 



ranging from 0° to 60° (A7 = 6°) is imposed in both 
calculations. 

The moments of inertia T K are calculated using the 
Inglis-Bclyaev formula [59, 60] 
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where k = 1,2,3 denotes the axis of rotation, and the 
summation of i, j runs over the proton and neutron quasi- 
particle states. The mass parameters B^ v (fi, 7) are given 
by [61] 
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The mass parameters B^ in Eq.(7) can be converted 
into the forms of Bpp, Bp^, S 77 by using the following 
relationships [53], 
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where the coefficients ao2 = aoo/v / ^22 = a oo/2, with 
a o = 9r%A w / 3 /WTt 2 , and r = 1.2. 



The potential V co \\ in Eq.(l) is given by, 

y coll (/3, 7 ) =^tot(/3,7)-AKib(/3,7)-AK-ot(/3,7), (10) 

where Etot is the total energy from constrained mean- 
field calculations. The Al^ib and Al/ ro t are the zcro- 
point-cncrgy of vibration and rotation respectively, 



AKib(/?,7) = \Tt[M^M {2) ], 
AK.ot(/?,7) = J2 AV W(A7), 



(11) 

(12) 



where the matrix M.t n ) is determined by Eq.(8) with in- 
dices fi 1 v running over and 2 and corresponding mass 
quadrupole operators defined as Q20 = 2z 2 — x 2 — y 2 and 
Q22 = x 2 — y 2 . Moreover, AV flfl (f3 7 j) is calculated by 

lM(2),^(/3,7) 



4^ (3 ),^(/?,7) ; 



(13) 



where M( n ) tllv (j3, 7) is also determined by Eq.(8), but 
with the intrinsic components of quadrupole operator de- 
fined as. 



!2n 



-2iyz, fx = 1 

— 2xZ, jJL = — 1 

2ixy, \x = —2 



(14) 



The details about the solution of constrained 
RMF+BCS and SHF+BCS equations have been given 
in Refs. [33] and [62] respectively. 
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III. RESULTS AND DISCUSSION 
A. Collective parameters in Hamiltonian 

Figure 1 displays the moment of inertia along x- 
direction, mass parameters and squared proton radii in 
5DCH for 96,98,ioo Sr and 98,ioo,i02 Zr isotopes as func _ 

tions of axial deformation parameter /3 determined from 
the mean-field calculations with both SLy4 and PC-PK1 
forces. We note that the collective parameters in 5DCH 
do not change much in 96,98,ioog r and 98, 100, 102^5- j so _ 

topes at all the /3-7 deformation regions, except for the 
Bpp with /3 ~ 0.4 and 7 = 60° as well as the B 11 with 
j3 ~ 0.5 and 7 = 0°. Moreover, the parameters in 5DCH 
determined from both SLy4 and PC-PK1 mean-field cal- 
culations are quite similar. For the mean squared proton 
radii r^, the SLy4 predicts slightly larger value than the 
PC-PK1 by - 0.5% systematically. 

B. Spectroscopy of nuclear low- lying states 

Figure 2 displays the electric quadrupole transition 
strength B(E2 : 2f -> Of), excitation energy E x (2f) of 
2f state and the ratio i? 4 / 2 (= E x (4f )/E x (2f)) as func- 
tions of neutron number in Sr and Zr isotopes from the 
5DCH calculations with both the SLy4 and the PC-PK1 
forces. Systematically the B(E2 : 2f -> Of) (or E x (2f )) 
values from the 5DCH calculations with both forces are 
increasing (or decreasing) with the neutron number up 
to TV = 62 in both Sr and Zr isotopes, which indicates 
the existence of transition from spherical to prolate de- 
formed shapes. Compared with the experimental data, 
the evolution of B(E2 : 2\ — > Of ) and R4/2 with respect 
to the neutron number around -/V = 60 is much more 
dramatic (moderate) in the results of SLy4 (PC-PK1) 
calculations. Similar evolution behavior as that by the 
PC-PK1 force is also observed in the 5DCH calculations 
with the Gogny force D1S [52]. However, this dramatic 
decreasing in E x {2~i) from N = 58 to N = 60 in Sr 
and Zr isotopes, together with the sudden increasing of 
E x (2f) in 96 Zr is not reproduced in both calculations. 

C. Deformation energy surfaces and collective 
wave functions 

To understand the evolution character of spectroscopic 
quantities around TV = 60 shown in Fig. 2, we plot the 
deformation energy surfaces of 96,98,ioog r and 98,ioo,i022 r 

isotopes in plane from the constrained mean-field cal- 
culations with both SLy4 and PC-PK1 forces in Figs. 3 
and 4 respectively. In both calculations, there is always 
a weakly oblate deformed minimum coexisting with a 
prolate minimum in the deformation energy surfaces of 
96,98,ioo Sr and 98,100,102^ Furt hermorc, the absolute 

minimum in both calculations is shifted from the oblate 
to the prolate side as the neutron number increases from 



N = 58 to N — 60. However, the subtle balance be- 
tween these two minima is quite different in these two 
calculations, as shown in Fig. 5, where the total en- 
ergy in 96,98,ioo Sr and 98,ioo,i02 Zr as a function f axial 

deformation parameter /3 is plotted. The evolution of 
deformation energy curves, i.e., from weakly oblate de- 
formed 96 Sr, 98 Zr to large prolate deformed 98 Sr, 100 Zr in 
SLy4 calculations, is much more rapid than that in PC- 
PK1 calculations. Moreover, the oblate and prolate min- 
ima in the deformation energy curves of 96 > 98 > 100 Sr and 
98,ioo,i022 r are q U jte close in energy. A very large mixing 
of oblate and prolate configurations is expected in their 
ground states. In contrary to the PC-PK1 calculations, 
the SLy4 force gives a broader or deeper oblate minimum 
in 96 Sr and 98 Zr respectively and deeper or broader pro- 
late minimum in 98,ioog r and ioo,i022 r respectively. In 
other words, the 96 Sr and 98 Zr are more oblate, while 
98,ioo Sr and ioo,i02 Zr are mQre palate in the SLy4 cal- 
culations. 

The structural change in low-lying states is illustrated 
clearly in Fig. 6, where the distribution of the squared 
collective wave functions pj a in the f3-j plane for the 
Of and 2f states from the 5DCH calculations with both 
the SLy4 and the PC-PK1 forces is plotted. The pi a is 
defined as 

pi a (P, 7) = E \KkW, 7)| 2 /5 3 | sin 3 7 |, (15) 

K 

which follows the normalization condition, 

/ m\ d 1Pla ((3 ll ) = l. (16) 
Jo Jo 

Here, ^ K (f3, 7) is the collective wave function that cor- 
responds to the solution of 5DCH. The distributions of 
the squared collective wave functions in these two calcu- 
lations are quite different. The SLy4 predicts a sudden 
transition from oblate 96 Sr to very good prolate 98,ioog r ^ 
while, the PC-PK1 gives a coexistence picture for the 
ground states of 96,98,ioog r ^ w jth the dominant compo- 
nent changing from oblate to prolate and back to oblate 
again moderately. For the 2f state in the calculations 
with the SLy4 force, the dominant component in 9 Sr 
is oblate and it becomes well prolate in 98 Sr. While 
in the calculations with the PC-PK1 force, the domi- 
nant component is already prolate in 96 Sr. This pic- 
ture provides a interpretation for the evolution char- 
acter of B(E2 : 2+ -> Of), E x (2+) and R i/2 shown 
in Fig. 2. Compared with the experimental data for 
B(E2 : 2+ -> 0+) and R i/2 , the SLy4 (or PC-PK1) force 
overestimates (or underestimates) somewhat the slop of 
shape transition around N = 60. 

D. Isotope shifts and monopole transition 
strengths 

Nuclear charge radii or isotopic shifts are good indica- 
tors of shape changes along isotopic chains. In Rcfs. [32] 
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FIG. 2: (Color online) The electric quadrupole transition strength B(E2 : 2f — > 0^), excitation energy ^(2^") and the ratio 
i? 4 / 2 (= E x (4f)/E x (2^)) as functions of neutron number in Sr and Zr isotopes from the 5DCH calculations with (left panel) the 
SLy4 force and (right panel) the PC-PK1 force, in comparison with the experimental data [63]. The R4/2 values for vibration 
(2.00) and rotation (3.33) limits are indicated with horizontal dotted lines. 
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FIG. 3: (Color online)The deformation energy surfaces of 

96,98,100g r and 98,100,102 Zr isotopes m the ^ planei from the 

constrained SHF+BCS calculations with the SLy4 force for 
the ph channel and density-dependent 5 force for the pp chan- 
nel. All energies are normalized to the absolute minimum. 
Each contour line is separated by 0.5 MeV. 



and [33], the evolution of charge radii with the neutron 
number in Sr and Zr isotopes around N ~ 60 has been 
studied with the self-consistent constrained mean-field 
calculations. The charge radii corresponding to differ- 
ent local minima in the deformation energy surface are 
compared with the data. It has been shown that a rapid 
change in nuclear shape is essential to reproduce the ex- 
perimental charge radii. However, in these two studies, 
the beyond mean-field correlation effect on nuclear charge 
radii has not been examined. In Rcf. [64], the dynamic 



FIG. 4: (Color online) Same as Fig. 3, but from the con- 
strained RMF+BCS calculations with the PC-PK1 force for 
the ph channel and a separable force for the pp channel. 



quadrupole correlation effect on charge radii of a large set 
of even-even nuclei has been studied in the framework of 
configuration mixing of projected axially deformed states 
based on a topological Gaussian overlap approximation. 
It has been shown that the dynamical correlation leads 
to an overall increase of radii, and it might also reduce 
the charge radii for some specific nuclei. Therefore, it is 
interesting to show the evolution of charge radii with the 
correlation effect from quadrupole fluctuation. 

Figure 7 displays the isotope shifts corresponding to 
different configurations in Sr and Zr isotopes from the 
mean-field calculations with both SLy4 and PC-PK1 
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FIG. 5: (Color online) The total energy in 96 . 98 > 100 Sr and 
98,100,102^ ag a f unc tion D f ax i a i deformation parameter f3 
from the constrained mean-field calculations with both the 
SLy4 and the PC-PKI forces. All energies are normalized to 
the value at = 0. 
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FIG. 6: (Color online) The distribution of squared collective 
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96,98,ioo Sr from the 5DCH C ai cu i a tions with both the SLy4 
and PC-PKI forces. 
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FIG. 7: (Color online) The isotope shift (normalized to 
N = 50) in Sr and Zr isotopes as a function of neutron number 
from the mean-field calculations of both SHF+BCS with the 
SLy4 force and RMF+BCS with the PC-PKI force. The cor- 
responding 5DCH calculated results and experimental data 
from Refs. [65, 66] are given as well. The "obi.", "pro.", 
"sph." are used to label the character of mean-field configu- 
rations. 



forces. In addition, wc present the 5DCH predicted iso- 
tope shifts, in which, the quadrupole fluctuation effect 
with triaxiality is included. It is shown clearly that the 
overall of isotope shifts could be reproduced much bet- 
ter by the 5DCH calculations with dynamic correlation 
effect. The shape evolution picture consistent with that 
exhibited in the B(E2 : 2^ — > Of), i? 4 / 2 as well as the 
distribution of squared wave functions is shown again in 
the isotope shifts, a sudden rising of which from N = 58 
to N = 60 due to the transition from weakly oblate de- 
formed to large prolate deformed shapes is slightly over- 
estimated (or underestimated) in the 5DCH calculations 
with the SLy4 (or PC-PKI) force. 

Electric monopole transition strengths are a model- 
independent signature of the mixing of configurations 
with different mean-square charge radii [69] and there- 
fore provide a good way to illustrate the occurrence of 
shape-coexistence. Figure 8 shows the evolution of elec- 
tric monopole transition strength p 2 (E0 : Of — > 0^), 
given by the off-diagonal clement of charge radius, 



p 2 (E0;0+ ^0+) = 



(0 2 + |E fc e fc r£|0+ 



(17) 



with Rq = 1.2 A 1 ' 3 fm, with respect to the neutron 
number in neutron-rich Sr and Zr isotopes. The sys- 
tematic of p 2 (E0 : Of —> Of) is similar in the calcula- 
tions with the SLy4 and PC-PKI forces. Both calcula- 
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FIG. 8: (Color online) The electric monopole transition 
strength p 2 (E0 : — > 0+) from the 5DCH calculations with 
both the SLy4 and PC-PK1 forces, in comparison with the 
experimental data [67, 68]. 



tions predict the same peak position, i.e., at 96 Sr and 
100 Zr, which are the nuclei before and after the dramatic 
transition respectively. Quantitatively, however, as ex- 
pected from the distribution of squared wave functions 
for 96,98,ioo Sr in Fig _ g the SLy4 ( or PC-PK1) force 

predicts a much weaker (or stronger) mixing of oblate 
and prolate configurations in their ground states, and 
therefore smaller (larger) p 2 (E0 : 0^ — > Of ) values in 
96 ' 98 ' 100 Sr. Compared with the experiment data for 98 Sr, 
the SLy4 (or PC-PK1) force underestimates (or overes- 
timates) the p 2 (E0 : O2 — > 0^) value by a factor of two 
approximately. Similar evolution behavior is also shown 
in Zr isotopes. However, except for 98 Zr, the difference 
in magnitude between the results of SLy4 and PC-PK1 
calculations is smaller than that in Sr isotopes. 



E. Single-particle energy levels and spin-orbit 
splitting 

The nuclear low-lying states reflect information about 
the underlying single-particle structure and vice versa. 
Therefore, to understand the evolution character of col- 
lectivity in Sr and Zr isotopes around N = 60 in the 
calculations with both the SLy4 and PC-PK1 forces, in 
Fig. 9, we plot the neutron and proton single-particle en- 
ergy levels in 96 < 98 Sr as functions of the axial deformation 
parameter (3. In general, the minima in the deforma- 
tion energy curve are associated with a shell effect due 
to the low level density around the Fermi energy. It is 
shown in Fig. 9 that, for the proton single-particle en- 



ergy levels, both calculations predict a large shell gaps 
around the Fermi energy level in the prolate side, which 
is mainly formed by two levels split from the degener- 
ate 7rl(7 g/ /2 state due to the deformation (or Jahn- Teller) 
effect. The size of this energy gap does not change too 
much from 96 Sr to 98 Sr in both calculations. It means 
that the proton plays a minor role in the rapid nuclear 
shape transition from the mean-field point of view. For 
the neutron single- particle energy levels, both calcula- 
tions predict two evident shell gaps around the Fermi 
energy level in oblate and prolate sides, however, the de- 
tails of the single-particle structure are quite different. 

Compared with the PC-PK1 force, the SLy4 force pre- 
dicts a stronger spin-orbit splitting for neutrons (by a 
factor of ~ 1.1) for all the states, as shown in Fig. 10, 
where the splitting of neutron spin-orbit doublet states, 



21 + 1 



j k =£±l/2, 



(18) 



in the spherical states of 96 Sr and 98 Sr as a function 
of the orbital angular momentum £ is plotted. e n ij < is 
the energy of single-particle state with quantum num- 
bers (n, £,j<). Consequently, the position of vlgr/2 state 
is pushed up and the vlhn/2 state is pulled down com- 
pared with the PC-PK1 results as shown in Fig. 9. As 
a result, the shell gaps around the Fermi energy at the 
spherical point and the minima of the deformation en- 
ergy curves are quite different. For the spherical point 
(/3 = 0), a relatively large shell gap at N = 56, which 
provides a mechanism responsible for the observed much 
higher E x (2f) in 96 Zr, is shown in the SLy4 calculation, 
but not in the PC-PK1 calculations. Moreover, com- 
pared with the PC-PK1 results for 96 Sr, the position of 
v2d§/2 state is almost the same, but the ulhn/2 state is 
much lower in the SLy4 results. As a result, the shell gap 
around the Fermi level in the prolate side, mainly formed 
by the K = 3/2 component of v2d§/2 orbit, the intruded 
K = 1/2 component of v2f 7 / 2 orbit, and other two lev- 
els with K = 3/2,5/2 split from fl/in/2 orbit, is much 
smaller than that in the PC-PK1 calculations. In 98 Sr, 
the energy of vlhn/2 state is shifted up, which broadens 
the shell gap significantly in the prolate side. This big 
change in the energy gap around Fermi energy is respon- 
sible for the sudden onset of large prolate deformation 
at N — 60 given by the SLy4 calculations. In the PC- 
PK1 calculations, however, this change in the shell gap of 
prolate side is more moderate. We note that, similar as 
the PC-PK1 results, the shift of ulh 11/2 state in 96 < 98 Sr 
from the Gogny D1S calculations is small and the change 
in the shell gap of prolate side is not significant [70]. 
On the other hand, compared with the SLy4 force, the 
PC-PK1 predicts a larger shell gaps in the oblate side 
of 96 ' 98 Sr, which provides the mechanism responsible for 
the strong mixing of prolate and oblate shapes in their 
ground states. 
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FIG. 10: (Color online) Splitting of neutron spin-orbit doublet 
states [c.f.(18)] in the spherical states of 96 Sr (open symbols) 
and 98 Sr (filled symbols) as a function of the orbital angular 
momentum £ from the mean-field calculations with both the 
SLy4 (triangle) and the PC-PK1 (circle) forces. 



IV. SUMMARY 

In summary, the rapid structural change in low-lying 
collective excitation states of neutron-rich Sr and Zr iso- 
topes has been studied by solving a 5DCH with pa- 
rameters determined from both the RMF and SHF cal- 
culations. Pair correlations are treated in the BCS 
method with either a separable pairing force or a density- 
dependent zero-range force. The isotope shifts, excita- 
tion energies, electric monopolc and quadrupole transi- 
tion strengths have been calculated and compared with 



corresponding experimental data. The calculated re- 
sults with both the PC-PK1 and the SLy4 forces ex- 
hibit a picture of spherical-oblate-prolate shape transi- 
tion in neutron-rich Sr and Zr isotopes. However, com- 
pared with the experimental data, the PC-PK1 (or SLy4) 
force predicts a more moderate (or dramatic) change in 
most of the collective properties around N = 60 and a 
much stronger (or weaker) mixing between oblate and 
prolate configurations in their ground states. The dif- 
ference between these two calculations is mainly because 
of the quite different structure in neutron single-particle 
states, mostly caused by the different spin-orbit inter- 
action strengths. Moreover, the sudden broadening of 
neutron shell gap in the prolate side, mainly formed 
by the K = 3/2 component of ^2d 5 / 2 , the intruded 
K =1/2 component of ^2/7/2, and other two compo- 
nents of ulhn/2 state, has been shown to be responsible 
for the rapid shape transition at N = 60. However, it has 
to be pointed out that the rapid change in the excitation 
energy of the first 2 + state has not been reproduced in the 
calculations with both the PC-PK1 and SLy4 forces. In 
particular, even though the SHF+BCS calculation with 
the SLy4 force indeed predicts a sizable spherical shell 
gap at N = 56, the corresponding 5DCH calculation is 
not able to reproduce the suddenly increased excitation 
energy for the first 2 + state at 96 Zr. A further beyond 
mean-field investigation is required. 
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